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У статті досліджуються проблеми та існуючі методи оптимізації фінансових показників мережевих 
аптек. Запропоновано генетичний метод з модифікацією оператора мутації для вирішення даної проблеми. 
Принципова відмінність розробленого генетичного алгоритму, порівняно з існуючими аналогами, полягає в 
можливості керування математичним розподілом значень рішення, що дозволяє запобігти передчасній збіж- 
ності генетичного алгоритму та використати всі запропоновані гени у долях згідно з обраною моделлю розпо- 
ділу. Крім цього, запропоновано три модифікації генетичного оператора мутацій. Перша модифікація полягає 
у виборі значень замін для мутації не випадковим чином, а з ряду, що підкоряється закону нормального розпо- 
ділу. Друга модифікація служить для визначення доцільності мутації хромосоми, спираючись на знання ретро- 
спективних та прогнозних даних з використанням прогнозної моделі штучної нейронної мережі (ШНМ). Третя 
модифікація полягає у комбінації двох зазначених модифікацій. Застосування розроблених методів призведе 
до більш ефективного використання площі аптек, зменшення незадоволеного попиту та, в кінцевому резуль- 
таті, до зменшення роздрібної вартості ліків за рахунок зменшення видатків на зберігання та обслуговування 
неоптимально завантажених площ аптеки. 

Ключові слова: генетичний алгоритм, еволюційний алгоритм, оптимізація фінансових показників, 
мінімізація часу перебування товару на складі 
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Вступ 

У сучасному конкурентному середо- 
вищі швидкість та правильність прийнят- 
тя рішень є ключовим фактором успіш- 
ності рітейлера, яким 1 є аптечна мережа. 
Головними індикаторами успішності ро- 
боти аптеки є її фінансові показники: при- 
буток та оборот. 


Впливати на ці показники можна різ- 
ними методами, але одним з найефектив- 
ніших є оптимізація асортименту складу 
аптеки. 

Питання оптимального асортименту 
є важливим як для аптеки, що працює три- 
валий час, так і для аптеки, що тільки буде 
відкриватись. 
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На вибір асортименту впливають такі 
фактори як площа аптеки, собівартість ме- 
дикаментів, оборотність асортиментних 
позицій, маркетингові фактори, фактори 
сезонності, екологічні фактори, фактори 
географічного положення аптеки відносно 
місцевої інфраструктури 1 т. ін. 

Оптимізація асортименту призведе 
до більш ефективного використання пло- 
щі аптек, зменшення незадоволеного по- 
питу та, в кінцевому результаті, до змен- 
шення роздрібної вартості ліків за рахунок 
зменшення видатків на зберігання та об- 
слуговування неоптимально завантажених 
площ аптеки. 

У ході роботи було розроблено гене- 
тичний метод багатокритеріальної оптимі- 
зації з модифікацією оператора мутації 
для дослідження ступеня впливу факторів 
на фінансові показники аптек та вибору 
моделі оптимізації. У роботі досліджу- 
ються проблеми та існуючі методи опти- 
мізації фінансових показників мережевих 
аптек, розробляється генетичний метод з 
модифікацією оператора мутації для вирі- 
шення проблеми управління асортимен- 
том аптечної продукції, заснований на 
еволюційних методах. У результаті вико- 
ристання запропонованих методів і засо- 
бів, виявляються найбільш значущі, З точ- 
ки зору підвищення ефективності, проце- 
си прийняття рішень у системі управління 
асортиментом. 

Постановка завдання 

Прибуток аптеки залежіть від таких 
базових факторів: (1) сума продажу; 
(2) валовий прибуток; (3) довжина складу; 
(4) показники захворюваності населення. 

Математично це можна виразити як 
функцію від багатьох змінних (1): 


воув/Об Х» Ха Ха). а) 


де К(х) - умовний фінансовий показник; 
хі - значення суми продажів; х2 -значення 
валового прибутку; хз - значення довжини 
складу аптеки; х/ - Значення показників 
здоров'я від обсягів викидів забруднюю- 
чих речовин у повітря. 

У цій роботі синтезується модель за- 
лежності фінансових показників від впли- 


ву асортименту продукції аптечної мере- 
жі. Перевіряється гіпотеза щодо впливу 
процентного співвідношення між групами 
товарів в асортименті на прибутковість 
аптеки. 

Аналіз літературних даних та 
постановка проблеми 

Оптимізація асортименту в аптеці - 
це комплекс заходів, спрямований на кіль- 
кісну і якісну зміну Його структури з ме- 
тою підвищення раціональності і ефектив- 
ності. Основна мета - гармонізувати асор- 
тимент, що призведе до мінімізації витрат 1 
збільшення прибутку аптеки. 

Поширеним способом оптимізації 
асортименту є аналіз дефектури, що дозво- 
ляє виявити попит на відсутні препарати 
шляхом фіксації попиту покупців перво- 
стольником в обліковій системі або на па- 
перових носіях (журнал незадоволеного 
попиту або обліку дефектури). Проте, він 
має і свої недоліки, оскільки можливе спо- 
творення даних за рахунок недостатньо 
ефективної системи відстежування дефек- 
тури. Окрім цього, не завжди можна отри- 
мати об'єктивні дані, оскільки дослідження 
споживчих переваг тільки в одній торговій 
точці часто не є репрезентативним | 11. 

У роботі |2| запропоновано АВС- 
аналіз для проведення повноцінного аналі- 
зу асортименту. Досить часто в аптеках для 
дослідження асортименту використовують 
АВС-аналіз. Ідея методу АВС-аналізу бу- 
дується на підставі принципу Парето: «за 
більшість можливих результатів відповідає 
відносно невелике число причин», зараз 
більш відомого як «правило - 20 на 30». 
Даний метод аналізу отримав значний 
розвиток завдяки своїй універсальності 1 
ефективності. 

За допомогою цього аналізу товари 
розбиваються за ступенем впливу на за- 
гальний результат. Причому принципом 
угруповання може бути величина вируч- 
ки, що отримується від конкретної групи 
продуктів, обсяг продажів або будь-які ін- 
ші параметри. Часто виручка більш пока- 
зова як критерій угруповання. Угрупован- 
ня за обсягом продажів може бути адек- 
ватним у тому випадку, якщо аналізовані 
товари однорідні за складом 1 ціною. 
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Робота |3| надає опис використаного 
ХУУ7-аналізу, який допомагає оцінювати і 
порівнювати стабільність продажів товар- 
них груп або окремих товарів різного типу 
попиту або різних цінових категорій. За- 
стосовується для оптимізації товарних 
запасів 1 визначення частоти замовлення 
товару. 

ХУУаналіз використовує показник, 
що характеризує потреби аптеки в товар- 
ному запасі. 

Алгоритм проведення ХУ/-аналізу 
включає наступні етапи: 

1) вибір об'єктів аналізу (товарна група, то- 
варна одиниця, постачальники, клієнти і 
т.п.); 

2) визначення параметрів аналізу (одиниці 
продажів, обсяг продажів, дохід, серед- 
ній запас, кількість замовлень 1 т.п.) і пе- 
ріоду аналізу (тиждень, місяць, квартал, 
півріччя, рік); 

3) визначення коефіцієнтів варіації для 
аналізованих ресурсів; 

4) угруповання ресурсів відповідно до 
зростання коефіцієнта варіації; 

5) розподіл за категоріями Х, У, /; 

6) графічне - представлення результатів 
аналізу. 

Необхідно представити задачу так, 
щоб її розв'язання можна було б записати 
у вигляді генотипу, тобто вектору значень 
(генів). Оптимальною буде така стратегія 
управління асортиментом, яка мінімізує 
суму всіх витрат, пов'язаних зі створен- 
ням, зберіганням і нестачею запасів, в 
одиницю часу або за певний (у тому числі 
і нескінченний) проміжок часу. У найза- 
гальнішому вигляді завдання управління 
асортиментом зводиться до знаходження 
такого розміру асортименту хі у момент 
часу І, який мінімізує загальну функцію 
витрат (2): 


Бі) -» тіп. 2) 


Управління асортиментом полягає у 
відшукуванні такої стратегії поповнення 1 
витрати запасами, при якій функція витрат 
набуває мінімального значення. Проста 
модель управління запасами представлена 
таким чином. Нехай функції 4А(/), В(ї) ви- 
ражають, відповідно, поповнення запасів, 
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їх витрату за проміжок часу |0, 1). Рівень 
запасу у момент г визначається основним 
рівнянням запасів (3): 


кі) ЕУНА)-ВІО), 0) 


де Ко-- початковий запас у момент 1-0 |7. 

Аналіз (1-3| за предметною областю 
дозволяє стверджувати, що розробка ме- 
тодів оптимізації фінансових показників є 
доволі актуальною задачею. Досліджено 
проблеми та існуючі методи управління 
асортиментом мережевих аптек. На основі 
досліджених методів аналізу асортименту 
можна зробити висновок про те, що комп- 
лексні методи аналізу асортименту пока- 
зують необхідність розгляду цілої групи 
показників ефективності асортименту. 
Склад цих показників і їх вплив на під- 
сумкову оцінку відрізняється залежно від 
особливостей асортименту, самої аптеки 1 
сформованої ринкової кон'юнктури. Тому 
методи |2-5| повинні не тільки вибирати, 
але й адаптуватися, виходячи з їх поточної 
ситуації на конкретній аптеці. Проаналізу- 
вавши методи |6-7| управління асортимен- 
том, можна зробити висновок, що лінійні 
моделі не мають повного набору можли- 
востей для вибору варіантів структури 
асортименту, оскільки вони дозволяють 
отримати оптимальне рішення тільки на 
один плановий період і не розглядають 
його зв'язок із показниками попереднього 
і наступного періоду. Тому для оптиміза- 
ції фінансових показників мережевих ап- 
тек та підвищення ефективності їх діяль- 
ності доцільно використовувати більш 
складні нелінійні, зокрема нейромережеві, 
моделі, які дозволять  апроксимувати 
складні багатовимірні нелінійні залежнос- 
ті з високою точністю. 

Мета і задачі дослідження 

Мета дослідження - створення гене- 
тичного методу багатокритеріальної опти- 
мізації з модифікацією оператора мутації 
для оптимізації фінансових показників 
мережевих аптек. 

Для вирішення завдання оптимізації 
фінансових показників аптек, таких як 
функції від розподілу товарних груп в 
асортименті аптеки, необхідно визначити 
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такий процентний розподіл груп товарів, 
при якому будуть виконуватись наступні 
умови: (1) максимізація прибутку аптеки; 
(2) мінімізація часу перебування товару на 
складі («довжини складу»); (3) покращен- 
ня прогнозованих показників прибутку 
аптеки; (4) варіабельність асортименту ап- 
теки, що забезпечує представленість мак- 
симальної кількості товарних груп в асор- 
тименті аптеки. 

Розробка модифікації оператора 
мутації при розв'язанні задачі оптиміза- 
ції фінансових показників аптек 

Задача оптимізації, що виникла і яка 
розв'язується в цій роботі, характеризу- 
ється великою кількістю змінних і, як на- 
слідок, великим обсягом простору пошу- 
ку, що не дає можливості дослідити все 
різноманіття рішень за прийнятний час. У 
зв'язку з цим виникла проблема практич- 
ної можливості розв'язання цієї задачі оп- 
тимізації: знайти ефективний або хоча б 
досить простий у практично важливих ви- 
падках алгоритм її вирішення. Для розв'я- 
зання цієї задачі було прийнято рішення 
використовувати еволюційні методи (ЗІ, 
які, порівняно з методами повного перебо- 
ру, дозволять скоротити обчислювальні 
витрати і розв'язати задачу оптимізації 
швидше та ефективніше. 

Генетичні алгоритми |9)| на цей мо- 
мент є найбільш відомим представником 
еволюційних методів оптимізації. 

Один із недоліків відомих еволюцій- 
них алгоритмів (10) полягає у відсутності 
механізму врахування обмежень оптиміза- 
ційного завдання. 

Тому, для усунення зазначеного не- 
доліку, у цій роботі запропоновано гене- 
тичний метод багатокритеріальної оптимі- 
зації з модифікацією оператора мутації 
для вирішення завдання оптимізації. Роз- 
глянемо деяку систему, що складається з 
двох підсистем, які описуються багато- 
екстремальними рівняннями: 


2угеУщ 4) 


25-25(Х2, У2 р (5) 


де ХІ, у, Х», у2 - параметри системи; 21, 


Сформулюємо проблему розв'язання 
задачі багатоцільової оптимізації такої 
системи: 

ж жо Ж 

РО ЗР(2|,22 /- тах Р(хУрХ2.У2), (6) 
де  Р- комплексна цільова функція; 
Хітіп З ХІ 2 Хітах) Хатіп 2 Х2 2 Хотах) 
Уїтіп 2 У 2 Уїтах) Удтіп З У2 2 Удтах- У зазначе- 
ному випадку функція Р є, по суті, компо- 
нентом багатоцільового показника якості 
Р(21,22) 1 перетворює сукупність таких 
компонентів | у | скалярний цільовий 
показник (11-12. 

Одним із найбільш поширених підхо- 
дів до врахування обмежень є метод штраф- 
них функцій |13|, основна ідея якого поля- 
гає в тому, що придатність індивіда обчис- 
люється не тільки залежно від відповідного 
йому значення цільової функції, а й від міри 
порушення обмежень: 


(7 
Лтез5 (Х -Дх) - б 'ЛМИ): 5 то 
7-1 


де і- номер покоління; 9-1, якщо вирішу- 
ється завдання мінімізації; 0--1, якщо вирі- 
шується завдання максимізації; /(х) -штраф 
за порушення )-го обмеження; /) -дійсне 
число /(1)-(С3)" (14-16). 

У вибраному методі штрафних функ- 
цій обчислення значення / (х) відбувається 
динамічно, залежно від міри порушення 
обмежень, за формулою для /-й ітерації, а 
значення /(1) -(С.1)": 


тах Зо, РАКУ, а (8) 


Л (х) - з 
| п (9 


,/егаіт 


де 2(х).0, п(х)-0 - обмеження задачі | 171. 
Отже, формула (7) придатності інди- 
віда має вигляд: 


(9) 
т 
Іпезз («Ла МС) - 5 7? (х.. 
7-1 

Перевагою методу динамічних штра- 

фів (13-20) є те, що він вимагає набагато 
менше параметрів, ніж інші методи штраф- 
них функцій. Замість вибору з набору фік- 
сованих рівнів порушення обмежень у 
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цьому методі штраф розраховується 
динамічно. 

Для покращення якості роботи алго- 
ритму та розширення його можливостей 
було розроблено модифіковані оператори 
мутації. 

Перша модифікація оператора мута- 
ції полягає у наступному: нові значення ге- 
нів для модифікації обираються не як ви- 
падкове число, а з ряду випадкових чисел, 
що підпорядковуються закону нормально- 
го розподілення |20-24). 

Алгоритм когнітивно-стильової 
детермінації виконується у наступній 
послідовності: (1) обрання хромосоми для 
мутації; (2) генерація умовно-випадкового 
масиву чисел за законом нормального роз- 
поділення, що за розміром приблизно до- 
рівнює розміру хромосоми. Медіанна точ- 
ка та середньоквадратичне відхилення роз- 
поділення співпадає з медіанною точкою 
ряду; (3) вибір генів для мутації випадко- 
вим чином; (4) заміна генів на значення з 
нормалізованого | випадкового ряду; 
(5) повернення хромосоми у популяцію. 

У хромосомі Азага2...ар вИПадковим 
чином вибирається | К позиція (біт) 
1«К п. Далі проводиться інверсія зна- 
чення гена в К позиції: ак-ак | 251. 

У когнітивно-стильовій детермінації 
значення гена після оператора мутації роз- 
раховується за формулами: 


ї с; - (б; -сідтприха0 (10) 
че с; - д(,Б, -а,) прих «1! 
Ь 
І І 
є 
б(уу)є уч  ТАХ/ ), 10 


де х - ціле випадкове число, що набуває 
значення 0 або 1; Є|0,1| - випадкове 
дійсне число; Єртах-- Максимальна кількість 
епох алгоритму; Б - параметр, що задається 
дослідником |261. 

Крім того, якщо протягом досить 
великого числа поколінь не відбувається 
збільшення пристосованості, то застосову- 
ються «мала» і «велика» мутації покоління. 
При «малій» мутації покоління до всіх 
осіб, крім 1095 кращих, застосовується опе- 
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ратор мутації. При «великій» мутації кож- 
на особа або мутує, або замінюється на 
таку, яка випадково згенерувала (27-29). 

Алгоритм Моейс (інтелектуальної) 
мутації полягає у використанні ШНМ у 
процесі мутації. Однією з цілей запропо- 
нованої модифікації є забезпечення тільки 
«позитивної» мутації, тобто такої, яка по- 
кращує фенотип хромосоми. Такий опе- 
ратор виконується у такій послідовності: 
(1) обрання | хромосоми для мутації; 
(2) застосування стандартної мутації; 
(3) використання раніше навченої ШНМ 
для прогнозування прибутковості, пере- 
давши на вхід мережі «оригінальної» хро- 
мосоми та «мутованої»; (4) порівняти від- 
повіді ШНМ і, тільки якщо «мутована» 
хромосома забезпечує кращу прибутко- 
вість, додати її в популяцію. В іншому 
випадку, в популяцію повертається 
«оригінальна» хромосома (301. 

Математично штучний нейрон, зазви- 
чай, представляють як деяку нелінійну 
функцію від єдиного аргументу - лінійної 
комбінації всіх вхідних сигналів. Цю функ- 
цію називають функцією активації або 
функцією спрацьовування, передавальною 
функцією |311. 

Функціонування | нейрона | можна 
описати формулою: 


М 

В и; ер 
М 

нео Ні «р 


12) 


де у - вихідний сигнал нейрона; у//...ууу - СИ- 
наптичні вагові коефіцієнти; и/...му -ВХІДНІ 
сигнали ШН; у - порогове значення (20). 

Модель (13) може бути представлена 
у вигляді: 


(13) 


де уоє У, йо-і. 

Третім типом модифікації оператора 
мутації є комбінація попередніх двох 
методів. 

Алгоритм Мегеег (об'єднавчої) мута- 
ції виконується у такій послідовності: 
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-обрання хромосоми для мутації; 

-застосування «нормалізуючої» мутації; 

-використання раніше навченої ШНМ для 
прогнозування прибутковості, передав- 
ши на вхід мережі «оригінальної» хро- 
мосоми та «мутованої»; 

-порівняти відповіді ШНМ 1, тільки якщо 
«мутована» хромосома забезпечує кра- 
щу прибутковість, додати її у популя- 
цію. В іншому випадку, в популяцію 
повертається "оригінальна" хромосома 
1321. 

Як функцію активації /, що наведена 

у формулі (13), можна використовувати 

функцію (14): 


Ііхеб п14) 


у на РО 


Окрім виразу (14), як функції актива- 
ції / також можуть застосовуватися порого- 
ві функції |22| вигляду (15) та (16) 


, Іі,хеб (15) 
о 
І, хі (16) 
Хо) «1-1х«-і 
х | єї 


Розроблено генетичний метод багато- 
критеріальної оптимізації з модифікацією 
оператора мутації, який переважає по на- 
дійності і швидкості порівняно з методами 
повного перебору. Окрім цього, модифіко- 
ваний генетичний алгоритм, наділений 
способами врахування обмежень, є ефек- 
тивним інструментом для розв'язання зада- 
чі оптимізації асортименту в аптеці. Це, у 
свою чергу, призведе до більш ефективно- 
го використання площі аптек, зменшення 
незадоволеного попиту та, в кінцевому ре- 
зультаті, до зменшення роздрібної вартості 
ліків за рахунок зменшення видатків на 
зберігання та обслуговування неоптималь- 
но завантажених площ аптеки. 

Експерименти методів ініціалізації 
початкової популяції еволюційного 
алгоритму 

Вибірка даних містила інформацію 


про структуру чека клієнтів аптечної мере- 

жі. Структура чека характеризується інфор- 

мацією про наявність певних товарів у ньо- 

му. Основними ознаками (атрибутами), що 

характеризують структуру чека покупця є: 

-х1 - ідентифікатор клієнта (КеуСизіотег) 
- унікальний номер, що дозволяє одно- 
значно визначити конкретного клієнта 
аптечної мережі; 

-Х2- 5КООТУ - кількість найменувань то- 
варів у чеку; 

-х3 - АТСОПу - кількість АТС груп; 

-х4 - 5аІс55ит - сумарна вартість товарів 
у чеку; 

-Х5 - МагвіпЗит - торгова точка, яка здійс- 
нює продаж товару; 

-Х6-- ОгдегОбу - кількість чеків певного 
клієнта на момент поточної покупки; 

-ху- АугОгдег - середній чек (грошова 
оцінка) певного клієнта; 

-Х8 - ОгдегКом8Оу - розмір знижки; 

-Хо-- АугРовійоп5ОГгу - середня кількість 
товарів у чеку певного клієнта; 

-Х10- - структура товарів у чеку, що пред- 
ставлена у вигляді: 


ЧІ б С5; " ; ат) 


де і - найменування і-го товару в чеку; 

Сі- кількість і-го товару в чеку; 5; -вар- 

тість одиниці /-го товару в чеку. 

Для вирішення завдання оптимізації 
фінансових показників аптек, як функції 
від розподілу товарних груп в асортименті 
аптеки, необхідно визначити такий відсот- 
ковий розподіл груп товарів, при якому 
будуть виконуватись наступні умови: 
-максимізація прибутку аптеки; 
-мінімізація часу перебування товару на 

складі («довжини складу»); 

-покращення прогнозованих показників 
прибутку аптеки; 

-варіабельність асортименту аптеки, що 
забезпечує представленість максималь- 
ної кількості товарних груп в асорти- 
менті аптеки. 

За одиницю групи товару було 
обрано групу класифікації довідника 
лікарських засобів  Моріон (Анатомо- 
терапевтично-хімічна класифікація - між- 
народна система класифікації лікарських 
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засобів. Використовуються скорочення: 
латиницею АТС (від Апаїотісаї 
Трегареціїс Сретісаї)). 

Важливою вимогою до даних для мо- 
делювання є їх якість. Якщо дані містять 
т.з. «шум», сезонну компоненту, викиди, 
пропуски - це негативно впливає на точ- 
ність прогнозів та якість моделей. 

При обробці пропусків даних порож- 
ні значення замінюються на медіанне зна- 
чення, розраховане за формулою (18): 


ІТ ма м 


Ме 


де Хме - нижнє значення медіанного інтер- 
валу; ім-- медіанний інтервал; 5ме- сума 
спостережень, що була накопичена до по- 
чатку медіанного інтервалу; /мг - КІЛЬКІСТЬ 
спостережень у медіанному інтервалі (26). 

Таким чином забезпечується міні- 
мальна статистична погрішність від зна- 
чень ряду. 

Для нормалізації даних використову- 
ється метод нормалізації МіпМах (23). 
Нормалізоване значення змінної х розрахо- 
вується за формулою (19) 


заія х - тіп (х) 
тах (х) - тіп (х)| (19) 


Одним із методів оцінки ефективності 
сформованої корзини асортименту є забез- 
печення прибутку у прогнозному періоді. 
Для надання такої оцінки доцільно викорис- 
товувати поширені методи прогнозування. 

Прогнозування можна здійснювати 
за абсолютним значенням або за напрям- 
ком тренда. 

До методів прогнозування за абсо- 
лютним значенням відносяться |291: 

- лінійна регресія; 
-поліноміальна регресія. 

Для прогнозування за напрямком, 
тренд можна представити як набір дискрет- 
них класів 30). У найпростішому випадку 
бінарної класифікації виконується умова 
(20). 
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0, Д/х) «0 


РО91 у да)»о" 


(20) 


де Р(х) - прогнозований клас; ЛДЙх) - зміна 
значення щодо попереднього періоду. 

Для оцінки якості прогнозних моде- 
лей використовується Інформаційний кри- 
терій Байеса (Вауевіап іпбоглайоп сгіїегіоп 
(ВІС)), середня абсолютна помилка (Меап 
АБяоїше Етгог (МАБ)) та середньоквадра- 
тична помилка (Меап 5ашаге Ептгог (М5Е)) 
132). 

Як навчальні моделі, використовува- 
лись дані про співвідношення груп АТС в 
асортименті аптеки. 

Шари ШНМ наведені в таблиці І. 


Таблиця І. Шари ШНМ для прогнозу- 
вання руху тренду 


Тип шару Кількість Функція 
нейтонтв активатті 
Вхідний (перший) 14 Кеї 0 
Прихований 23 Веї, 
(пуутий ) 

Прихований (третій) 20 Веї 0 
Прихований 10 Ве 0 
(четветтий) 

Вихідний (п'ятий) 1 Ве. 0 


У таблиці 2 наведені порівняльні 
характеристики моделей прогнозування, 
розраховані на даних продажів аптеки за 
період 2013-2017 років. 


Таблиця 2. Якісні характеристики прог- 
нозованих моделей 


Модель ВІС МАЕ | М5Е 
Лінійна - 102,8 1 - , 
пегпестя 
Поліноміальна | 6357 - | 
петтесія ' 

Логістична 94,99 | | 
пегпестя 
ШНМ - 0,024 | 0,1067 
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Складність оптимізації в цьому ви- 
падку полягає в тому, що оптимізувати 
треба не один параметр, а два (прибуток та 
оборотність товару), 1 при цьому забезпе- 
чити присутність усіх груп товарів в асор- 
тименті аптеки. Якщо прибуток виразити, 
як Р, а довжину складу через І, то у загаль- 
ному вигляді функція для оптимізації 
набуде вигляду |371: 


ПРІнЕ-, о 


де Е - фітнес функція. 

На рис. | представлений типовий 
розподіл АТС-груп товарів у розрізі дов- 
жини складу (оборотності товару). На 
цьому шррисунку зображена залежність 
середньомісячного прибутку на одиницю 
продукції від середньої довжини складу. 
Можна побачити, що за обсягом аптечних 
продажів лідує група АТС406, зокрема має 
такі показники, як середньомісячний 
прибуток, який складає 0.60 та середню 
довжину, що дорівнює 0.4. 


з | АТС406 
о 


АТСАДІІЗСО5 


АТС4232 АТС4238 
У АТС4234 


АТСА229 АТО 
ато ТС4035 


й прибуток 
1-1 


0.50 


Ус3280 
АТС3611 


АТСИ71 
АТС3892 
Т 


Середньомісячни 


Т Т Т 
03 04 0.5 06 07 
Середня довжина склада, нормалізована 
Рис. 1. Типовий розподіл АТС-груп товарів у 
розрізі довжини складу (оборотності товару) 


Можна припустити, що класичний 
ГА швидко зійдеться на очевидному рі- 
шенні розподілити більший відсоток за 
групами, що мають найбільшу прибутко- 
вість та невелике значення довжини скла- 
ду. На рис. 2 зображено відсотковий розпо- 
діл груп товарів класичного ГА. Як можна 
побачити, ГА розподілив 94 95 на групи 
АТС406 та АТС4232, а на інші групи 
всього 6 90. 
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- 


і АТС-група 
Рис. 2. Результат роботи класичного ГА 


Але для аптеки такий розподіл є не- 
прийнятним, тому що при цьому сформу- 
ється незадовільний запит покупців, що 
призведе до зниження кількості відвідува- 
чів та, як наслідок, зниження прибутку. 

Тому модифікований ГА повинен 
виконувати додаткову умову - розподіл 
груп товарів має бути близький до нор- 
мального розподілу (функція Гауса) |3581. 

Для вирішення цього завдання в ці- 
льову функцію було введено поняття 
штрафу алгоритму. Це означає, що якщо 
оцінювана популяція не відповідає нор- 
мальному розподілу, то значення фітнес- 
функції падає пропорційно. Також були 
введені вагові коефіцієнти для прибутку 1 
довжини складу, що дозволило гнучко 
керувати пріоритетом показників (391. 

Рівняння (21) набуло вигляду (22): 


2Рчей (22) 


І з 
е 


ГДОР,І) є 


де Р - очікуваний прибуток; І, - довжина 
складу; А - ступінь близькості розподілен- 
ня до нормального. 

Величина А - це результат тесту 
Андерсона Дарлінга, що показує, наскіль- 
ки розподіл вибірки співпадає з нормаль- 
ним розподілом. Чим більша ця величина, 
тим швидше зростає фітнес-функція. 

Після коригування фітнес-функції ГА 
показав наступні результати (рис. 3). На 
рис. 3 зображено відсотковий розподіл 
груп товарів класичного ГА з використан- 
ням штраф-алгоритму: як можна побачити, 
розподіл долі товарних груп в асортименті 
став більш варіабельним. 


О ЄМ. Федорченко, А.О. Олійник, О.О. Степаненко, С.К. Корнієнко, А.С. Харченко, Д.А. Гончаренко 65 


ло 


12596 


"ль 
| ле 
я 93 
82 
тб» 
о 
589 
51 

Зм ль 

305 

45 
09р- - 


АТС-група 
Рис. 3. Розподіл груп товарів 
(штраф-алгоритму) 
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На рис. 4 представлено результати ро- 
боти запропонованого методу когнітивно- 
стильової детермінації, а саме відсотковий 
розподіл груп товарів. Як можна побачити, 
розподіл долі товарних груп в асортименті 
став більш близьким до нормального розпо- 


ділу (функція Гауса). 
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АТС-група 
Рис. 4. Розподіл груп товарів 
(когнітивно-стильова детермінація) 


Наступна модифікація оператора му- 
тації полягає у внесенні в роботу оператора 
мутації «інтелекту». Ця модифікація слу- 
жить для визначення доцільності мутації 
хромосоми, спираючись на знання ретро- 
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спективних та прогнозованих даних з вико- 
ристанням ШНМ як прогнозовану модель. 
Результати роботи модифікованого опера- 
тора представлені на рис. 5. 

На рис. 5 зображено відсотковий 
розподіл груп товарів Моейс-мутації, на 
якому можна побачити, що, порівняно з 
когнітивно-стильовою детермінацією, роз- 
поділ долі товарних груп в асортименті 
Моепіс-мутації не підкорюється закону нор- 
мального розподілу. Але, при цьому, моди- 
фікований ГА забезпечує досить прийнятні 
параметри вихідної множини. 
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Рис. 5. Розподіл груп товарів 
(Моейс-мутація) 


Третя модифікація полягає у комбі- 
нації двох зазначених модифікацій. Резуль- 
тати роботи модифікованого оператора 
представлені на рис. 6. На рис. б зображено 
відсотковий розподіл товару при об'єд- 
нанні когнітивно-стильової детермінації та 
Моепіс-мутації, на якому видно, що розпо- 
діл долі товарних груп в асортименті став 
більш близьким до нормального розподілу 
порівняно з Моейс-мутацією, але менш 
близьким, ніж у  когнітивно-стильовій 
детермінації. 
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Андерсона-Дарлінга, середньоквадратична 

- ЕД похибка, очікуваний середній приріст при- 

и бутку, середня довжина складу нормалізо- 

ль вана. Під тестом Андерсона-Дарлінга |22| 

5 я вважається критерій, який призначений 

жу для перевірки, наскільки розподіл вибірки 

58» співпадає ї з нормальним розподілом. 

4 Середньоквадратична похибка є середньо- 

квадратичною відстанню між прогнозова- 

ним і фактичним значенням. Критерій при- 

бутковості розраховувався як різниця між 

"АННИ ТП - рівнем валового доходу торгівлі та рівнем 
З Зо З витрат обігу. 

У таблиці 3 наведено порівняльні да- 

ні роботи класичного ГА та модифікова- 

них версій. 
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Оцінка ефективності 0 генетичних 
алгоритмів проводилася за параметрами: 
тривалість розрахунку ітерацій, тест 


Таблиця 3. Порівняльний аналіз модифікацій ГА 


Ж Тривалість Тест Середньо- ная за Середня 
і Я Тип моделі розрахунку, | Андерсона - | квадратична ренане довжина склада 
ітерацій, С Дарлінга похибка при бутку, 9 нормалізована 
1 | Лінійна регресія 875 0,8878315 0,8080 0,50525923 0,4955 
з; 1 ЗНОНІШОМІНЛВНа 983 0,8971466 0,998 0,5212312 0,5397 
регресія 
9 ЗТОГОШНО 787 0,8772425 0,977 0.509569 0,46345 
регресія 
4 ШНМ 991 0,8878415 0,967 0,5345789 0,39445 
5 Класичний ГА 994 3,67510"7 0,9358 0,5217804 0,3964 
6 | ГА зі штрафом 974 0,9978415 0,754 0,5407766 0,44987 
Когнітивно- 
7 стильова 393 0,9837505 0,782 0,5564371 0,4587188 
детермінація 
8 Моегіс-мутація 804 0,9986964 0,638 0,5397302 0,4946176 
9 | Мегеег-мутація 1024 0,9976742 0,685 0,5438465 0,4544578 
Обговорення нейромоделі складає 1024с, при викорис- 


Можна побачити, що при викорис- 
танні Моейс- та Мегеег-мутацій збільшу- 
ється тривалість розрахунку (зокрема, при 
використанні Мегеег-мутації час синтезу 


танні Моейс-мутації -- 304с порівняно з ча- 
сом 737с при використанні логістичної ре- 
гресії). Таке збільшення зумовлене вико- 
ристанням ШНМ, яка тренована на реаль- 
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них співвідношеннях, що не підкорюються 
закону нормального розподілу. Але при 
цьому ГА забезпечує достатньо оптимальні 
параметри вихідної множини. 

Метод | «Мегеег-мутація» показує 
більш успішні результати порівняно з логіс- 
тичною регресією, а саме за показником 
тесту Андерсона-Дарлінга (при використан- 
ні Мегеег-мутації показник тесту складає 
0,9976742 на відміну від логістичної 
регресії, у якій показник тесту Андерсона- 
Дарлінга дорівнює 0,3772425) можна ствер- 
джувати, що розподіл вибірки є більш прий- 
нятним до нормального розподілу. Отже, 
фітнес функція  Мегеег-мутації зростає 
швидше порівняно з логістичною регре- 
сією. А також даний метод забезпечує кра- 
щу прибутковість (при Мегеег-мутації при- 
бутковість становить 0,543846596, а при 
логістичній регресії 0,50956995), проте 
вимагає більше обчислювальних та часових 
витрат (при використанні Мегеег-мутації 
час синтезу нейромоделі складає 1024с) з 
причини, покладеної в його основу як ідеї 
об'єднання когнітивно-стильової детермі- 
нації та Моейс-мутації. 

При порівнянні Моебіс-мутації з кла- 
сичним ГА можна зробити висновок, що 
Моейс краща за показниками тривалості 
розрахунку (зокрема, при використанні 
Моейс-мутації час синтезу нейромоделі 
складає 804с порівняно з часом 994с при 
використанні класичного алгоритму) та 
середнім приростом прибутку (показник 
прибутку | у  Моейс-мутації складає 
0,539730295, а у класичного ГА становить 
0,521780495), але поступається за показни- 
ками довжини складу (зокрема, при вико- 
ристанні Моепс-мутації складає 0,4946176, 
а у класичного ГА становить 03964). 
Порівнюючи Моегіс-мутацію з ГА зі штра- 
фом, що є також розвитком класичного 
ГА, можна сказати, що ті алгоритми 
рівноправні за ефективністю. 

Найкращі результати | показала 
когнітивно-стильова детермінація, тобто 
вона дає виграш у тривалості розрахунку, 
що становить 393с на середньому прирос- 
ті прибутку, який складає 0,556437190, що 
набагато більше своїх аналогів. 

Отримані результати дозволяють зро- 
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бити висновок про те, що запропонований 
підхід до вирішення завдань оптимізації 
фінансових показників мережевих аптек 
дає змогу збільшити середній приріст при- 
бутку аптеки на (0,539о та мінімізувати 
довжину складу до 0,45 з урахуванням 
середньоквадратичної похибки 0,685. 

Таким чином, у роботі запропонова- 
но 1 обгрунтовано новий підхід до вирішен- 
ня завдання оптимізації процесу роботи 
відділу закупівлі ліків, основна ідея якого 
полягає у використанні модифікованого 
генетичного методу для оптимізації пара- 
метрів моделі з контролем математичного 
розподілення значень вихідної хромосоми, 
з метою підвищення ефективності (стій- 
кості) ГА як еволюціонуючої системи. Роз- 
роблений алгоритм з ефективності вирі- 
шення задачі оптимізації на множині тесто- 
вих даних перевершує методи регресії та 
класичний ГА. 

Висновки 

Наукова новизна роботи полягає у 
тому, що запропоновано модифікований 
генетичний метод для оптимізації пара- 
метрів моделі з контролем математичного 
розподілення значень вихідної хромосо- 
ми. Принципова відмінність розробленого 
генетичного алгоритму, порівняно з існу- 
ючими аналогами, полягає в можливості 
керування математичним розподілом зна- 
чень рішення, що дозволяє запобігти 
передчасній збіжності генетичного алго- 
ритму та використати всі запропоновані 
гени у долях згідно з обраною моделлю 
розподілу. Крім цього, запропоновано три 
модифікації генетичного оператора мута- 
цій. Перша модифікація полягає у виборі 
значень замін для мутації не випадковим 
чином, а з ряду, що підкоряється закону 
нормального розподілу. Це дозволяє отри- 
мати більш «нормалізовані» хромосоми та 
прискорити збіжність генетичного алго- 
ритму з завданим законом розподілу. 
Друга модифікація служить для визначен- 
ня доцільності мутації хромосоми, спи- 
раючись на знання ретроспективних та 
прогнозованих даних з використанням 
ШНМ як прогнозованої моделі. Це дозво- 
ляє направляти мутацію, не додаючи в 
популяцію гени, які не покращать вихідну 
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популяцію. Третя модифікація полягає у 
комбінації двох зазначених модифікацій. 

Порівняльний аналіз роботи класич- 
ного ГА та модифікованих версій показує, 
що найкращі результати досягнуті у 
когнітивно-стильовій детермінації. Крім 
того, аналіз показав, що Моебйіс- та Мегеег- 
мутації істотно ефективніші за середнім 
приростом прибутку, незважаючи на вели- 
ку тривалість розрахунку. Отже, отримані 
результати дозволяють зробити висновок 
про те, що запропоновані модифікації опе- 
ратора мутації є доцільними та ефектив- 
ними для вирішення проблеми оптимізації 
асортименту аптеки. Застосування розроб- 
лених методів призведе до більш ефектив- 
ного використання площі аптек, зменшен- 
ня незадоволеного попиту та, в кінцевому 
результаті, до зменшення роздрібної вар- 
тості ліків за рахунок зменшення видатків 
на зберігання та обслуговування неопти- 
мально завантажених площ аптеки. 
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І. Кедогспепко, А. ОЇйіпук, 
А. 5'ерапепко, 5. КогпіїепКо, 
А. Крагсрепко, Р. СопспагепКо 

Деуеїортепі ої еуоішбііопагу аїієсо- 
гійт5 іп Фе ргобіетя ої тапа?іпе Ше 
а55огітепі ої пебуогК ррагтасіе5 

Тре агіїсіе ехріоге5 Ше ргобіетя апа 
ехізитя, пеїодйз8 ої орійпиліпе Ше Нпапсіа! 
регіогтпапсе ої пекмуогк рбагтасіез. Ргоро- 
5ед а бепейс тефодй ут Фе птоаїйсайоп ої 
Фе тикайоп орегаог (о 50Іуе Ше ргобіет ої 
тапавіпє Ше гапеєе ої ррагтпасецйса! рго- 
Фисі5 Ба5зед оп еуоїийопагу теШодйз. ТПе 
Гпатепіа Фійегепсе Бебуеєп Ше Чеусіо- 
ред єепейс аїсогіййт апа 15 ехіз іп сойп- 
(еграгі8 15 Фе абішу іо сопігої Фе птаШета- 
цсаї фізіібицоп ої Фе уашез ої Фе 501ибоп, 
муісб ргеуепі5 ргеташге сопуегеепсе ої Ше 
гепейс аїєогіййт апа и5ез5 ай ргорозед 
депе5 іп Насйоп5 ассогдїпеє (о Ше спобеп 
фізіібицоп поадсеі. Га ададїоп, (ргее поадій- 
сайопя (о Ше репейс орегаїог ої плиайоп5 
аге ргоро5ей. Те Нг5: піодійсаноп 15 0 
споо5е Ше 5иб5ійшиоп уаїез їог а ппикайоп 
пої їп а гапдогл у/ау, бий йот а 5егіе5 Шаї 
обеує Ше Іам/ ої погплаї дФі5сібийоп. Тре 
зесопа птодійсайоп 5егуе5 (0 Фегегтипе Фе 
арргоргіагепезз ої Ше спготоботе птикайоп, 
рБазед оп Кпом/едєе ої геїго5ресйує апа 
ргедїсйуе даїіа цз5іпе а5 а ргедїсйуе поасі ої 
УММ. Тре іга тодійсайоп 15 Ше сопібі- 
пайоп ої Фе їуо птіодійсайоп5 тепйопед. 
Тре арріїсайоп ої (Ше деуеіоред пефодз 
МШ Ісай (0 а тоге ейесйуе иц5е ої Фе рраг- 
тасу агеа, іо гедисе иптеї детапа апа, цій- 
таїеу, їо гедисе Ше геїаї! собі ої дгиє5 Бу 
гедисіпє Ше со5іє ої 5ї0гіпє апа 5егуісіпе, 
Фе зоборйта! Іоадіпе ої Фе рбагтасу. 
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